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Solid  Oxide  Fuel  Cell  (SOFC)  integrated  into  Micro  Gas  Turbine  (MGT)  is  a  multivariable  nonlinear  and 
strong  coupling  system.  To  enable  the  SOFC  and  MGT  hybrid  power  system  to  follow  the  load  profile  accu¬ 
rately,  this  paper  proposes  a  self-tuning  PID  decoupling  controller  based  on  a  modified  output-input 
feedback  (OIF)  Elman  neural  network  model  to  track  the  MGT  output  power  and  SOFC  output  power. 
During  the  modeling,  in  order  to  avoid  getting  into  a  local  minimum,  an  improved  particle  swarm  opti¬ 
mization  (PSO)  algorithm  is  employed  to  optimize  the  weights  of  the  OIF  Elman  neural  network.  Using 
the  modified  OIF  Elman  neural  network  identifier,  the  SOFC/MGT  hybrid  system  is  identified  on-line,  and 
the  parameters  of  the  PID  controller  are  tuned  automatically.  Furthermore,  the  corresponding  decou¬ 
pling  control  law  is  achieved  by  the  conventional  PID  control  algorithm.  The  validity  and  accuracy  of  the 
decoupling  controller  are  tested  by  simulations  in  MATLAB  environment.  The  simulation  results  verify 
that  the  proposed  control  strategy  can  achieve  favorable  control  performance  with  regard  to  various  load 
disturbances. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

With  increasing  global  energy  needs  and  environmental  con¬ 
cerns,  interest  in  new  energy  sources  is  growing.  Solid  Oxide  Fuel 
Cell  (SOFC)  integrated  in  Micro  Gas  Turbine  (MGT)  is  a  promising 
and  preferred  method  for  generating  electric  power  in  the  future, 
having  a  high  efficiency  and  low  pollution. 

Transients  in  the  load  may  have  a  significant  impact  on  the  life 
of  the  SOFC/MGT  hybrid  power  system.  One  of  the  reasons  is  that 
load  transients  often  involve  significant  peaks  in  power  relative  to 
the  steady-state  load.  Therefore,  the  primary  control  goal  for  the 
SOFC/MGT  hybrid  system  is  to  allow  the  power  plant  to  deliver  the 
desired  power  output. 

There  have  been  very  few  reports  on  controlling  the  power  of 
the  SOFC/MGT  hybrid  system  during  the  last  several  decades.  In 
Refs.  [1,2],  a  dynamic  model  was  presented  for  control  of  the  inte¬ 
grated  SOFC  and  turbine  hybrid  system  and  the  power  supplied 
by  the  SOFC  system  was  controlled  by  manipulating  the  fuel.  In 
Ref.  [3],  an  unstable  power  output  was  observed  due  to  the  fluc¬ 
tuation  of  gas  composition  in  the  fuel.  A  specially  designed  fuel 
controller  succeeded  not  only  in  allowing  the  hybrid  system  to  fol¬ 
low  a  step  change  of  power  demand,  but  also  maintained  the  system 
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power  output  stable.  However,  the  SOFC/MGT  hybrid  system  is  a 
multivariable  and  strong  coupling  nonlinear  system.  This  coupling 
affects  the  effectiveness  of  a  specific  loop  controller  on  the  corre¬ 
sponding  output,  and  in  some  case,  may  become  serious  and  cause 
many  difficulties  to  the  control  system  design.  How  the  power  sys¬ 
tem  is  decoupled  and  practical  controllers  are  designed  does  not 
feature  currently  in  the  open  literature. 

In  recent  years,  various  linearization  and  decoupling  methods 
have  been  presented  to  handle  this  problem.  The  usual  methods 
are  based  on  differential  geometry  method,  which  solve  a  group 
of  differential  equations  and  linearize  nonlinear  system  with  state 
feedback  [4-6].  However,  the  SOFC/MGT  hybrid  system  is  an  uncer¬ 
tain  nonlinear  system.  It  is  very  difficult  to  describe  such  a  system 
with  a  precise  mathematical  model. 

Neural  networks  are  considered  as  an  attractive  structure  to 
establish  the  mathematical  relationship  of  the  dynamic  system 
based  on  the  input-output  data.  As  a  kind  of  recurrent  neu¬ 
ral  network,  the  OIF  Elman  neural  network  is  more  ascendant 
than  the  static  neural  network  such  as  BP  and  RBF  neural  net¬ 
work  on  the  dynamic  characteristic.  Thereby,  it  is  now  widely 
used  in  the  areas  of  system  identification,  nonlinear  control  and 
prediction  control  [7-10].  Furthermore,  based  on  the  modified 
OIF  Elman  neural  network,  an  adaptive  proportional-integral- 
derivative  (PID)  decoupling  control  scheme  is  proposed  to  control 
the  SOFC  output  power  and  MGT  output  power.  The  tradi¬ 
tional  PID  controller  is  widely  applicable  in  many  fields  due  to 
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Nomenclature 

Pj  gas  pressure  at  state  j  (atm) 

Wj  gas  mole  flowrate  at  state  j  (mol  s-1 ) 

Wjj  mole  flowrate  of  gas  i  at  state  j  (mol  s_1 ) 

Xj  molar  fraction  of  gas  i 

Xjj  molar  fraction  of  gas  i  at  state  j 

r  moment  of  inertia  ( kg  m2 ) 

R  gas  constant  (8.314 J  mol-1  K-1 ) 

Van  anode  volume  (cm3) 

Vca  cathode  volume  (cm3) 

i  SOFC  stack  current  (A) 

Q3  enthalpy  of  gas  i  at  state  j  (J  mol-1 ) 

st  turbine  pressure  ratio 

rjt  turbine  efficiency 

hj  j  gas  enthalpy  at  state  j  (J  mol-1 ) 

Tj  gas  temperature  at  state  j  (K) 

N  SOFC  number 

Pac  stack  AC  output  power  (kW) 

vst  stack  output  voltage  (V) 

Pst  stack  output  power  (kW) 

Pmgt  MGT  power  (kW) 

Pt  turbine  power  (kW) 

r\ac  DC/AC  conversion  efficiency 

ridc  DC/DC  conversion  efficiency 

rjm  turbine  mechanical  efficiency 

W\  air  flowrate  (mol  s_1) 

Wq  fuel  flowrate  (mol  s-1 ) 

n  rotational  speed  (rpm) 

Tjr  compressor  isentropic  efficiency 

sr  compressor  pressure  ratio 

Ts  SOFC  stack  temperature  (I<) 

Pr2  consumed  power  by  the  fuel  compressor  (kW) 

Subscripts 
CFI4  methane 

CO  carbon  monoxide 

C02  carbon  dioxide 

N2  nitrogen 

02  oxygen 

H2  hydrogen 

FI20  water  (gas) 

i  state  i 

j  gas  species  j 


1 - r 


Anode  recirculation  SOFC 


Fig.  1.  Structure  diagram  of  SOFC/MGT  hybrid  power  system. 


sions  and  suggestions  for  future  work  are  finally  made  in  Section 
6. 

2.  Description  of  SOFC/MGT  hybrid  power  system 

The  heart  of  the  SOFC/MGT  hybrid  power  system  is  the  inter¬ 
nal  reforming  SOFC  (IR-SOFC)  stack,  which  is  designed  here  on  the 
basis  of  the  tubular  technology  pioneered  by  Siemens  Westing- 
house  in  the  1970s  [14].  In  order  to  avoid  excessive  thermal  stresses 
in  the  stack,  the  SOFC  operating  temperature  requires  preheating 
the  fuel  and  air  streams.  Thus,  it  is  mandatory  to  use  heat  exchang¬ 
ers  that  can  transfer  the  high  thermal  energy  rates  of  SOFC  outlet 
streams  increasing  the  air  and  fuel  inlet  temperature.  The  steam 
needed  to  support  the  internal  reforming  reaction  is  obtained  by 
re-circulating  part  of  the  anode  outlet  stream,  which  should  be 
controlled  to  supply  sufficient  steam  to  meet  the  steam  and  carbon 
ratio.  Therefore,  a  pressurized  SOFC/MGT  hybrid  power  system  is 
designed  in  this  paper,  shown  in  Fig.  1 .  It  consists  of  an  anode  recir¬ 
culation  SOFC  and  a  micro  gas  turbine.  The  operation  process  of  the 
SOFC/MGT  hybrid  system  has  been  described  in  the  Ref.  [25]. 

3.  Power  mathematical  model  of  SOFC/MGT  hybrid  power 
system 

Certain  assumptions  are  made  in  developing  the  dynamic  model 
for  the  SOFC/MGT  hybrid  generation  system.  These  assumptions  are 
all  valid  from  an  analytical  modeling  point  of  view. 


its  simplicity  and  robustness  [11-13].  However  the  SOFC/MGT 
hybrid  system  is  a  complex  coupling  nonlinear  system,  which 
can  hardly  be  solved  by  the  conventional  PID  controller  with 
fixed  parameters.  Thus  here,  the  self-learning  ability  of  the  OIF 
Elman  neural  network  is  used  to  tune  the  PID  parameters  on-line. 
The  coupling  is  treated  as  an  exterior  disturbance,  so  the  pro¬ 
posed  controller  is  used  to  eliminate  disturbance  and  improve  the 
tracking  performance.  A  power  mathematical  model  of  a  220  kW 
SOFC/MGT  hybrid  system  is  used  to  generate  the  data  required 
for  the  training  and  prediction  of  the  OIF  Elman  neural  network 
model. 

This  paper  is  outlined  as  follows:  Section  2  briefly  describes  the 
SOFC/MGT  hybrid  power  system.  In  Section  3,  a  power  mathemat¬ 
ical  model  of  the  hybrid  system  is  built  with  MATLAB  to  imitate 
the  real  220  kW  SOFC/MGT  system.  An  adaptive  PID  decoupling 
control  strategy  based  on  the  modified  OIF  Elman  neural  net¬ 
work  is  presented  for  power  decoupling  control  of  the  SOFC/MGT 
hybrid  system  in  Section  4.  The  responses  of  the  SOFC/MGT  hybrid 
system  to  certain  load  changes  are  studied  in  Section  5.  Conclu- 


•  Gas  mixtures  are  resolved  for  CH4,  CO,  C02,  H2,  H20,  N2,  and  02. 
All  other  species  are  assumed  negligible  for  thermodynamics. 

•  Both  air  and  fuel  are  assumed  as  the  ideal  gas,  which  makes  them 
satisfy  the  ideal  gas  state  equation. 

•  There  is  surplus  oxygen  available  for  complete  combustion  due 
to  the  fact  that  air  mass  flowrate  is  much  larger  than  the  fuel  mass 
flowrate  in  the  burner. 

•  No  heat  transfer  to  the  environment.  The  system  is  assumed  to 
be  well  insulated  from  the  environment. 

•  Each  cell  in  the  stack  is  assumed  to  operate  identically,  so  that 
a  single  SOFC  simulation  is  taken  as  representative  and  used  to 
calculate  full  stack  performance. 

•  Control  volumes  are  characterized  by  a  single  lumped  tempera¬ 
ture,  pressure,  and  species  mole  fractions  condition. 

3.1.  Valve  and  mixer  model 

In  order  to  supply  steam  for  the  reforming  process,  a  certain  part 

of  the  anode  exhaust  gas  must  be  recycled.  The  valve  is  used  for  the 


X.-J.  Wu  et  al.  /  Journal  of  Power  Sources  196  (201 1)  1295-1302 


1297 


control  of  fuel  and  gas  flows.  It  is  adiabatic  and  fluid  inertia  is  not 
considered.  Flowrate  and  molar  composition  of  the  input  and  the 
output  are  set  equal  at  any  load.  In  the  mixer,  the  amount  of  steam 
supplied  to  the  pre-reformer  is  determined  by  the  steam  carbon 
ratio  Xsc,  defined  as  the  ratio  between  the  supplied  steam  to  the 
cell  and  the  supplied  methane. 

3.2.  Pre -reformer  model 


The  mass  balance  for  the  SOFC  is  described  as  follows: 


Pi 0 Kin  ^10 ,i 
—RTS  dT 


—  W9  i  -  W10  j  +  Rr  i  +  Ran  j 


P5  Kca  ^5 ,j 

~RTh~dF 

where 


W4j-W5J  +  RcaJ 


(10) 

(11) 


Rr  =  [  —T\  r\  -  r2  r2  3ri  +  r2  -ri  -  r2  ] 


(12) 


The  supplied  fuel  is  reformed  to  hydrogen  and  carbon  monoxide  pan  =  [o  0  0  -r3  r3 
by  fuel  reforming  reaction  (Eq.  (1 ))  and  water-gas  shifting  reaction 
(Eq.  (2))  in  the  pre-reformer.  Rca  =  [  0  -0.5 r3  ] 


(13) 

(14) 


CH4  +  H20  CO  +  3H2 

(1) 

CO  +  h2o  co2  +  h2 

(2) 

The  reaction  rate  of  hydrogen  corresponding  to  the  electro¬ 
chemical  reaction  r3  is  directly  related  to  the  current  i,  which  is 
given  in  Eq.  (15): 


Therefore,  the  energy  balance  for  the  pre-reformer  is  given  by: 

Cre  W=  ^2{Ws  '  XiM  '  Ss’i)  ~J2{W9'  *9.''  ■  S9.‘)  +  Or  +  &  (3) 

i  i 

where  Oe  is  the  heat  capacity  of  reformer,  Qr  and  Qs  the  consumed 
heat  by  the  reforming  reaction  Eq.  (1)  and  the  shifting  reaction  Eq. 
(2),  which  can  be  calculated  as  follows: 

Qr  =  Pre,  1  *  (^9,CO  +  3frg,H2  -  h8,H20  ~  ^8,CH4)  (4) 

Qs  =  Pre, 2  •  (^9,C02  +  ^9,H2  “  ^8,H20  “  ^8,Co)  (5) 

Under  the  ideal  gas  supposition,  the  partial  molar  enthalpy  hjj 
is  calculated  by  the  formulation: 

ki  =  k 0+  [J  Cj/T)dT  (6) 

href 

where  j  =  l,  2,  •••,  16,  ig{CH4,  CO,  C02,  H2,  H20,  02,  N2}.  huo  is 
the  standard  enthalpy  of  gas  2,  and  Tref  is  the  reference  temperature 
(=298  K).  The  specific  heat  Cjti  of  gas  i  can  be  obtained  from  Ref.  [15]. 

3.3.  IR-SOFC  model 


Taking  into  account  ohmic,  concentration  and  activation  losses, 
applying  Nernst’s  equation  and  Ohm’s  law,  the  stack  output  voltage 
can  be  obtained  according  to  Refs.  [16-18]. 

The  output  power  of  the  SOFC  stack  Pst  can  be  obtained  by  the 
equation: 

Pst  =  ivst  (16) 

The  power  conditioning  unit  is  the  interface  between  the  SOFC 
output  DC  voltage  and  the  AC  grid,  which  consists  of  a  DC/DC 
converter  and  a  DC/AC  inverter.  Assuming  the  DC/DC  and  DC/AC 
conversion  efficiencies  are  pdc  and  ijac  respectively,  the  SOFC  output 
power  is  given  by: 

Pac  =  Pst  •  Pjdc  '  P}ac  (17) 

3.4.  Burner  model 

In  the  catalytic  combustion  burner,  the  residual  fuels  from  the 
SOFC  are  burned  away  which  increases  the  temperature  of  the  gas. 
The  following  reactions  are  being  considered  during  the  combus¬ 
tion: 


There  are  three  chemical  reactions  considered  in  the  SOFC  stack. 
These  reactions  are  the  methane  reformation  reaction  (Eq.  (1)), 
water-gas  shift  reaction  (Eq.  (2)),  and  electrochemical  reaction  (Eqs. 
(7)  and  (8)). 

H2 +  02“  ^  H20  +  2e“  (7) 

j02  +  2e~  02~  (8) 

Therefore,  to  calculate  the  SOFC  stack  temperature,  three  heat 
sources  should  be  determined,  namely:  (i)  heat  consumed  by  the 
reforming  reaction  Qi ;  (ii)  heat  consumed  by  the  shifting  reaction 
Q2;  (iii)  heat  generated  due  to  the  electrochemical  reaction  in  the 
SOFC  stack  Qs.  Assuming  that  temperature  is  uniform  in  a  stack,  the 
energy  balance  equation  for  the  SOFC  can  be  expressed  as: 


2H2+02  -*  2H20  (18) 

2C0+02  -*  2C02  (19) 

The  temperature  of  the  product  gas  can  be  determined  from: 

Cb  hf  =  ^ ^ + Wl2 '  ■  hu’i] + & 

i  j 

-J2w^xn,rhn,i  (20) 

i 

where,  i  g{CO,  C02,  H2,  FI20,  02,  N 2},j  g{H2,  CO},  Q,  is  the  combus¬ 
tion  heat  of  H2  and  CO  respectively. 

3.5.  Heat  exchanger  model 


Cs  ■  w  = W4  ■  ■  ^j) + W9  •  55*9.1  •  h,i) 

j  i 

-  W5  ■  5>5J  ■  h5J)  -  Ww  ■  ■  hio ,i)  -  Pst  O) 

j  i 

+  5J2*  O) 

k 

where  k  =  1,  2,  3,  ie{CH4,  CO,  C02,  H2,  H20},  je{N2,  02},  Cs  is  the 
stack  heat  capacity. 


Compact  heat  exchangers  are  often  used  in  recuperated  gas 
turbine  cycles.  As  one  kind  of  compact  heat  exchanger,  plate-fin 
heat  exchanger  is  chosen  in  this  model.  The  dynamics  of  the  heat 
exchanger  are  described  in  Ref.  [19]. 

3.6.  Compressor  model 

The  compressor  map  is  based  on  a  DLR  centrifugal  compressor 
map,  while  the  turbine  map  is  for  a  NASA-CR-1 74646  axial  turbine 
[20].  The  original  and  smoothened  data  files  of  these  turbomachines 
are  obtained  from  the  manufacturers.  According  to  the  similarity 
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Anode  exhaust 
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Inlet  gas  Outlet  gas 


Turbine  mdel 


Stack  current 

SOFC  power 

Anode  inlet 

Cathode  outlet 
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IR-SOFC  stack  model 

SOFC  power 


Fuel  compressor  power 


T urbine  power  Rotor  speed 


Burner  model 


d> 

Inlet  air 


Rotor  speed 

Power 

Inlet  air 

Outlet  air 

Air  compressor 


Outlet  cool  gas 


►  Inlet  hot  gas 


Outlet  hot  gas 


Inlet  cool  gas  Outlet  cool  gas 


Air  compressor  power 


Heat  exchanger  2  model 


MGT  power 


Heat  exchanger  1  model 
Fig.  2.  SOFC/MGT  hybrid  power  system  model  in  SIMULINK. 


theory,  the  features  of  the  compressor  can  be  described  as  a  func¬ 
tion  of  the  pressure  ratio,  reduced  mass  flow,  reduced  speed  and 
efficiency. 

When  the  operational  point  is  located,  pressure  ratio,  isentropic 
efficiency  and  flowrate  are  known,  the  exit  temperature  T2  of  the 
air  compressor  can  be  expressed  as  follows: 


3.8.  Shaft  model 


A  shaft  model  accounts  for  the  dynamics  of  the  rotating  mass  in 
the  gas  turbine  system,  which  is  modeled  as  follows: 


dn  900  n 

-J7  =  - n  •  l  MGT 

dt  nxiz 1 


(25) 


T2  =  7"!  +  —  ■  (4r,_1)/r'  -  1)  (21) 

hr 

where  adiabatic  coefficient  r\  can  be  obtained  from  Ref.  [21  ]. 

The  consumed  thermal  power  by  the  air  compressor,  Pr i ,  can  be 
calculated  from  the  following  equation: 

Pri  =  W\  •  (/i2  —  h\)  (22) 


3.7.  Turbine  model 

The  compressor  and  gas  turbine  are  assumed  to  be  coupled  on 
a  single  shaft.  As  a  consequence,  they  have  the  same  rotor  speed. 
The  maps  are  used  to  calculate  the  pressure  ratio  and  isentropic 
efficiency  as  a  function  of  mass  flowrate  and  rotor  speed  .The  ideal 
temperature  of  the  working  fluid  at  the  outlet  of  the  turbine  can  be 
evaluated  using  Eq.  (23). 

Ts  =  r14-7-14(i  (23) 

The  turbine  work  can  be  calculated  from  the  real  enthalpy 
change  and  the  flowrate,  which  is  expressed  as  follows: 

Pt  =  Wu-(hu-h15)  (24) 


where  r  is  the  moment  of  inertia.  Neglecting  the  generator  mechan¬ 
ical  loss,  the  output  mechanical  power  from  the  gas  turbine  PMgt  is 
defined  in  the  following  equation: 

PmGT  =  hm-Pt  -  Pn  -  Pr2  (26) 

The  power  mathematical  model  replaces  the  real  SOFC/MGT 
hybrid  power  system  to  generate  the  simulation  data  required  for 
the  identification  of  the  modified  OIF  Elman  neural  network  model. 
For  the  SOFC/MGT  dynamic  physical  model,  the  SOFC  output  power 
and  MGT  output  power  are  the  controlled  variables,  fuel  flowrate 
and  air  flowrate  are  chosen  as  manipulated  variables,  and  the  load 
current  is  considered  as  a  disturbance.  The  SOFC/MGT  hybrid  power 
system  model  developed  in  MATLAB  is  shown  in  Fig.  2.  The  param¬ 
eters  of  the  SOFC/MGT  hybrid  system  are  given  in  Table  1. 

4.  Power  decoupling  control  of  SOFC/MGT  hybrid  power 
system 

4.  \ .  Power  control  strategy 

The  structure  for  the  power  decoupling  control  of  the  SOFC/MGT 
hybrid  system  is  shown  in  Fig.  3,  which  combines  the  modified  OIF 
Elman  neural  network  and  the  PID  controller. 

In  Fig.  3,  m  =  1,  2,  rm(/<)  is  the  reference  power  (SOFC  out¬ 
put  power  and  MGT  output  power),  ym(k)  the  real  output  power, 
ymm(k)  the  modified  OIF  Elman  neural  network  identification 
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power,  em(k)  the  error  between  the  given  value  rm(k)  and  output 
ym(k )  in  every  instant,  um{k)  the  manipulated  variable.  According 
to  the  error  em(/<),  the  modified  OIF  Elman  neural  network  identi¬ 
fication  model  is  used  to  tune  the  parameters  of  the  conventional 
PID  controller  to  keep  the  system  stable.  For  two  inputs  and  two 
outputs  system,  the  coupling  effect  from  the  second  loop  is  treated 
as  exterior  disturbance  to  the  first  main  loop.  At  the  same  time, 
the  coupling  effect  from  the  first  loop  is  treated  as  exterior  distur¬ 
bance  to  the  second  main  loop.  Thereby,  the  PID  controller  based 
on  the  modified  OIF  Elman  neural  network  is  used  to  eliminate 
disturbance  and  improve  the  tracking  performance. 

4.2.  Modified  OIF  Elman  neural  network  model 

The  following  NARX  model  is  used  to  describe  the  controlled 
power  system: 

y{k)  =f[y{k  -  1), ...,  y{k  -  ny),  u{k  -  1), ...,  u(k  -  nu)]  (27) 

where  y  is  the  SOFC  output  power  and  MGT  output  power,  u  is 
the  fuel  flowrate,  air  flowrate  and  current,  ny  and  nu  are  the  lags 
of  the  output  and  input  respectively,  and  /( • )  is  a  nonlinear  func¬ 
tion.  In  this  section,  we  adopt  a  modified  OIF  Elman  neural  network 
identification  model  to  identify  the  nonlinear  function /( • ). 

The  OIF  Elman  neural  network  is  a  type  of  recurrent  neu¬ 
ral  network,  which  has  two  particular  layers  called  the  context 
layer  and  the  context  2  layer  besides  the  conventional  input, 
hidden  and  output  layers.  The  context  layer  and  context  2  lay¬ 
ers  are  used  to  memorize  the  former  values  of  the  hidden  and 
output  layer  nodes  respectively.  The  feed-forward  connections 
are  modifiable,  whereas  the  recurrent  connections  are  fixed. 
The  structure  of  OIF  Elman  neural  network  is  shown  in  Fig.  4 
[22]. 

In  Fig.  4,  wu  is  the  weight  between  the  input  layer  and  hidden 
layer,  wy  is  the  weight  between  the  hidden  layer  and  the  output 


Table  1 

Setting  parameters  used  in  the  model  of  SOFC/MGT  hybrid  system. 


Parameter 

Unit 

Value 

Cell  number  N 

- 

1152 

Active  area  A 

cm2 

834 

DC/DC  conversion  efficiency  r]dc 

- 

97% 

DC/AC  conversion  efficiency  r]ac 

- 

97% 

Compressor  pressure  ratio  er 

- 

2.9 

Compressor  isentropic  efficiency  r]r 

- 

78% 

Turbine  isentropic  efficiency  r]t 

- 

82% 

Steam/Carbon  ratio  Xsc 

- 

2.5 

Turbine  mechanical  efficiency  r]m 

- 

94% 

Moment  of  inertia  r 

kg  m2 

0.00127 

Stack  heat  capacity  Cs 

JK-1 

471 

Fig.  3.  Power  decoupling  controller  structure  of  the  SOFC/MGT  system. 


Fig.  4.  Structure  of  OIF  Elman  neural  network. 
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layer,  wc  is  the  weight  from  the  context  layer  to  the  hidden  layer, 
and  wyc  is  the  weight  from  the  context  2  layer  to  the  hidden  layer. 
xc{k)  and  x{k)  are  the  outputs  of  the  context  unit  and  the  hidden 
unit  respectively.  yc{k)  and  ym(/<)  are  the  outputs  of  the  context  2 
layer  and  output  layer  respectively,  a  and  y  are  the  feedback  gains 
of  the  self-connections  of  context  and  context  2  layers  respectively, 
0  <  a  <  1 ,  0  <  y  <  1 . 

The  mathematical  model  of  the  OIF  Elman  neural  network  is 


described  as  follows: 

x{k)  =f(wcxc{k)  +  wuu(k  -  1)  +  wycyc(/<))  (28) 

xc{k)  =  a  •  Xc(k  —  1 )  +  x(k  —  1 )  (29) 

yc{k)  =  y  ■  yc{k  - 1)  +  y(k  -  1)  (30) 

ym(k)  =  g(wyx(k))  (31) 

where  g(x)  is  often  taken  as  a  linear  function,  and  j[x)  is  often  taken 
as  a  sigmoid  function,  namely: 

/m=ttW  (32) 


The  standard  Elman  network  usually  adopts  BP  algorithm  to 
train  the  network’s  weights,  however  it  is  easy  to  get  locked  into 
local  minima  [23].  Therefore,  here  an  improved  particle  swarm 
optimization  is  adopted  to  train  the  weights  of  the  OIF  Elman  neural 
network  [24].  The  main  steps  can  refer  to  the  optimization  process 
of  Ref.  [25]. 


4.3.  Self-tuning  PID  decoupling  control  based  on  modified  OIF 
Elman  model 

The  digital  incremental  PID  control  algorithm  is  adopted  here, 


which  is  expressed  as  follows: 

tim{k)  =  um(k  —  1)  +  kpmxm(I)  +  kimxm(2)  +  kdmxm(  3)  (33) 

where 

)  =  0m(k)  —  em(k  —  1 )  (34) 

Xm(  2)  =  em{k)  (35) 

Xm{  3)  =  0m{k)  —  2em(k  —  1 )  +  em(k  —  2)  (36) 

where  kpm,  kim  and  kdm  are  the  proportional  factor,  integral  factor 
and  differential  factor  respectively. 

Define  the  performance  index  Hm  as  follows: 

Hm(fc)  =  \\emm2  =  T[rm(lc) -ym(k)]2  (37) 


Adopting  the  gradient  descent  method,  kpm,  kim  and  kdm  are 
adjusted  as  per  the  following: 

kPm(k )  =  kpm(k  —  1 )  —  tjp  — —  =  kpm  ( k  —  1 )  +  r/pem(k )-^- — xm(l ) 
Ol<pm  OUm 

(38) 


kim(k)  =  kim(k  kjm(k  -  1 )  +  /?,em(k)^xm(2) 

(39) 


kdm(k )  =  kdm(k  - 1 )  -  =  kdm(k  -  1 )  +  %em(k)^xm( 3) 

(40) 

where  rjp ,  rp  and  rjd  are  the  proportion,  integral  and  differential 
learning  rate  respectively,  9ym/9um  is  the  Jacobian  information  of 
controlled  object,  which  can  be  obtained  from  the  above  modified 
OIF  Elman  identification  results,  i.e.  dymldum  ^3yrnm/3um. 
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(b)  Response  of  the  SOFC  output  power  and  MGT  output  power 
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(c)  Response  of  the  SOFC  output  power  in  [100  Is,  1 101s] 


Fig.  5.  Response  of  the  SOFC/MGT  output  power  during  +15%  step  change. 


5.  Simulation 

In  this  section,  we  present  numerical  experiments  to  validate 
the  proposed  adaptive  PID  decoupling  control  scheme  based  on  the 
modified  OIF  Elman  identification  model  of  the  SOFC/MGT  hybrid 
system. 

From  the  above  mentioned  Sections  2  and  3,  firstly,  based  on  the 
mass  balance  and  enthalpy  balance  equations,  the  power  mathe¬ 
matical  model  of  the  SOFC/MGT  hybrid  system  described  in  Section 
2  is  built  with  MATLAB  to  imitate  the  real  220  kW  SOFC/MGT  hybrid 
system  in  Section  3.  With  the  mathematical  model,  the  training  data 
for  the  modified  OIF  Elman  neural  network  model  is  generated.  In 
order  to  obtain  available  identification  data,  the  input  signals  of 
the  power  physical  model  are  uniformly  random,  including  the  fuel 
flowrate,  air  flowrate  and  the  current.  A  set  of  600  data  is  collected 
from  the  simulation.  The  first  400  data  are  used  for  the  identifica¬ 
tion  of  the  modified  OIF  Elman  neural  network  identification  model, 
while  the  remaining  200  data  are  used  for  validation  purposes. 

In  this  study,  Kolmogorov  theorem  is  used  to  determine  the 
number  of  hidden  node,  and  is  given  as  2r+l,  where  r  is  num- 
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Fig.  6.  Response  of  the  SOFC/MGT  output  power  during  -15%  step  change. 


ber  of  input  units.  There  are  three  inputs  (air  flowrate  W\y  fuel 
flowrate  W6  and  current  i).  So,  the  hidden  layer  of  the  modified 
OIF  Elman  neural  network  is  calculated  to  be  7  nodes.  Here,  the 
improved  PSO  algorithm  is  applied  to  obtain  the  weights  of  the 
OIF  Elman  neural  network.  Therefore,  each  particle  consists  of 
3x7  +  7x2  +  7x7  +  2x7  parameters,  and  all  weights  are  chosen 
on  the  interval  [-1, 1].  The  parameters  of  the  improved  PSO  algo¬ 
rithm  are  set  as  follows:  swarm  size  m  is  40,  maximum  number  of 
iterations  Gmax  is  200,  accelerating  factors  C\  and  c2  are  1.8,  maxi¬ 
mum  inertia  weight  TxTmax  is  0.9,  minimum  inertia  weight  Tnm[n  is 
0.1,  step  size  of  the  inertia  weight  Atjj  is  0.05.  Using  MATLAB  to 
run  the  program  of  the  modified  OIF  Elman  model  and  after  200 
iterations,  the  optimal  weights  of  OIF  Elman  can  be  determined. 

Based  on  the  modified  OIF  Elman  neural  network  model,  the 
adaptive  PID  decoupling  controller  is  used  to  control  the  power 
of  the  SOFC/MGT  hybrid  system.  Firstly  assume  that  the  external 
load  of  the  SOFC/MGT  hybrid  system  has  a  +15%  step  change  at 
1 000  s.  The  simulation  time  is  3000  s,  and  the  sampling  time  is  0.1  s. 


The  proportion,  integral  and  differential  learning  rate  are  r\v  =  0.09, 
77i  =  0.3  and  ijd  =  0.2  respectively.  The  PID  parameters  kp,  /q  and  kd 
are  adjusted  by  self-learning  modified  OIF  Elman  neural  network 
until  the  error  approach  zero.  The  simulation  results  are  shown  in 
Fig.  5.  The  output  power  of  the  SOFC/MGT  hybrid  system  tracts  and 
decouples  the  reference  power  satisfactorily,  which  is  shown  in 
Fig.  5(a).  The  maximum  error  of  the  output  power  is  1 .02%,  and  the 
mean  steady-state  error  is  0.03%.  The  response  of  the  SOFC  output 
power  and  the  MGT  output  power  are  shown  in  Fig.  5(b).  The  spe¬ 
cific  change  of  the  SOFC  from  1001  s  to  1101  s  is  shown  in  Fig.  5(c). 
In  order  to  compare  with  the  traditional  PID  control,  now  assum¬ 
ing  the  external  load  of  the  SOFC/MGT  hybrid  system  has  a  -15% 
step  change  at  20  s,  the  simulation  results  are  shown  in  Fig.  6.  From 
Fig.  6,  the  traditional  PID  controller  is  incapable  of  controlling  the 
output  power  to  the  target  value.  However,  using  the  improved 
PID  decoupling  control  method,  the  SOFC/MGT  output  power  can 
effectively  tract  the  reference  power.  The  maximum  error  of  the 
output  power  is  1.58%,  and  the  mean  steady-state  error  is  0.06%.  In 
addition,  as  the  SOFC  output  power  is  the  core  output  power  of  the 
hybrid  system,  the  variation  curve  of  the  SOFC/MGT  output  power 
is  similar  to  the  SOFC  output  power. 

6.  Conclusions 

In  order  to  allow  the  SOFC/MGT  hybrid  power  system  to  deliver 
the  desired  power  output,  this  paper  employs  an  adaptive  PID 
decoupling  controller  based  on  a  modified  OIF  Elman  neural  net¬ 
work  identification  model  to  control  the  SOFC  output  power  and 
MGT  output  power.  The  data  required  for  the  training  and  pre¬ 
diction  of  the  modified  OIF  Elman  neural  network  identification 
model  is  generated  from  a  power  mathematical  model  of  a  220  kW 
SOFC/MGT  hybrid  power  system.  The  simulation  results  show  that 
the  proposed  adaptive  PID  decoupling  control  algorithm  based 
on  the  modified  OIF  Elman  neural  network  model  is  an  efficient 
method  to  solve  the  nonlinear  coupling  power  system  for  the 
SOFC/MGT  hybrid  system.  In  order  to  develop  an  integrated  control 
strategy,  in  the  future  other  transient  impacts  besides  current  will 
be  further  discussed  (e.g.,  ambient  temperature,  pressure  fluctua¬ 
tions  and  fuel  composition  changes). 
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